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Abstract. It has long been known that random regular graphs are with high probability good ex- 
panders. This was first established in the 1980s by BoUobas by directly calculating the probability that 
a set of vertices has small expansion and then applying the union bound. 

In this paper we improve on this analysis by relying on a simple high-level observation: if a graph 
contains a set of vertices with small expansion then it must also contain such a set of vertices that is 
locally optimal, that is, a set whose expansion cannot be made smaller by exchanging a vertex from 
the set with one from the set's complement. We show that the probability that a set of vertices satisfies 
this additional property is significantly smaller. Thus, after again applying the union bound, we obtain 
improved lower bounds on the expansion of random zi-regular graphs for Zi > 4. In fact, the gains from 
this analysis increase as A grows, a fact we explain by extending our technique to general A. Thus, in 
the end we obtain an improvement not only for some small special cases but on the general asymptotic 
bound on the expansion of Z\-regular graphs given by Bollobas. 
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1 Introduction 



A graph is called an expander graph if for all partitions of its set of vertices into two sets the ratio of the 
number of edges with endpoints in both sets over the size of the smaller set is bounded from below by 
a constant independent of the graph size. Roughly speaking, an expander is a very well-connected graph 
and we are usually interested in expanders which are regular and sparse. Expander graphs have long been 
an object of intense study in discrete mathematics and theoretical computer science. Over time, they have 
been used in a wide number of applications, such as for example error-correcting codes [16) . computational 
complexity [9115] . inapproximability for bounded-occurrence CSPs [14 and others too numerous to mention 
here (see [TU]). 

One of the most fundamental results in this area, established by BoUobas in 1988 [5], is that random 
Z\-regular graphs are expanders with high probability. Thus, though the existence of some expander graph 
families had already been known |6ll2j . Bollobas managed to show that in fact almost all Z\-regular graphs 
are expanders using surprisingly elementary methods (straightforward counting arguments, as opposed to 
the linear algebra techniques commonly used when dealing with expander graphs). Furthermore, Bollobas 
gave an asymptotic lower bound oi ^ — ©(v^) on the expansion of random Z\-regular graphs, and specific 
bounds for various small values of A. Interestingly, the asymptotic bound was later found to be optimal up 
to the coefficient of VA, as Alon showed that any Z\-regular graph (not just random graphs) must have a 
bi-section with expansion at most ^ — &{VA) [1]. 

The main objective of this paper is to revisit this classical result of Bollobas and to improve it through a 
more refined analysis. With this new analysis we will show higher lower bounds on the expansion of random 
Z\-regular graphs. 

Our results move in two main directions. First, we will obtain better lower bounds on the expansion 
of Z\-regular graphs for various small values of A. One of the reasons we are interested in this is because 
improvements in the currently known bounds could have applications in other areas, and in particular one 
example is the inapproximability of bounded occurrence CSPs. For instance, the fact that 6-regular graphs 
have expansion at least 1 (established in [5]) is invoked a number of times in this area (see |4I3I2| ). It is 
still open whether Z\ = 6 is the smallest degree for which this is true, and resolving this question would be 
of immediate use in inapproximability reductions. Of course, beyond this, finding what is the expansion of 
random Z\-regular graphs for small values of A is in itself a fundamental mathematical problem, intimately 
connected to a number of other basic graph-theoretic questions such as the size of the minimum bisection 
of random regular graphs. Unfortunately, little progress has been made since [5] and in most cases the 
bounds given by Bollobas are still the best known lower bounds, although some progress has been made in 
establishing corresponding upper bounds through results on the minimum bisection problem [131718] . 

Second, we will use our refined analysis to also improve the asymptotic bound on the expansion as a 
function of A, or more precisely, to show that the current bound is not tight. As mentioned, it is known that 
the correct value is ^ — 0{\/A) but the coefficients of y/A given by Bollobas in the lower bound (Vln 2 « 0.83) 
and Alon in the upper bound (j^^ ~ 0.13) are relatively far apart. Again, this is a fundamental question 
given as an open problem in [1 . But furthermore, the main reason we would like to give such a general 
analysis is to show that the refinements we make in the analysis of small values of A are more than just an 
ad-hoc trick that improves some small special cases, but in fact our technique applies to all A. 

To explain the high-level idea of our method let us first recall the proof of the main result of [5^ . The 
idea there is to calculate the probability that a certain fixed set of vertices of a random Z\-regular graph 
has exactly c edges connecting it to its complement. Clearly, as c becomes smaller that probability also 
becomes smaller, so the trick here is to calculate the largest value of c that still makes the probability 
o(2^"). Afterwards, an application of the union bound ensures that with high probability no set has c or 
fewer edges crossing the cut, giving a bound on the expansion. The best bound on c is found by directly 
writing out the probability that exactly c edges cross the cut and performing a tedious but straightforward 
calculation. 

Since the probability that c edges cross the cut is calculated exactly in [5 one might expect that the best 
place to look for an improvement may be in the application of the trivial union bound, which completely 
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ignores a very basic fact: The expansion ratios of two different sets of vertices are far from independent. 
Somewhat surprisingly, we give an argument that attempts to exploit this crucial fact while still relying on 
the union bound. The main idea is that if a set with small expansion exists, then there must exist a set with 
the same or smaller expansion which is locally optimal. A set is locally optimal in this context if exchanging 
a vertex from it with a vertex from its complement cannot decrease its expansion. Thus, to establish that a 
graph's expansion is above a certain value it suffices to show that no locally optimal set exists with expansion 
below this value. Again, we want to find the maximum c such that the probability that a fixed set of vertices 
is locally optimal and has c edges leaving it is o(2~") and, as might be expected, this probability turns out 
to be significantly smaller than the probability that the set simply has c edges leaving it, allowing us to end 
up with a larger c. 

The main technical challenge arising here is to calculate this probability accurately. In the case of [51 it 
is a simple counting exercise to determine the number of configurations where exactly c edges cross the cut 
obtaining a clean and relatively simple formula. In our case however, as we will argue, the local optimality 
condition imposes some constraints on the degrees of the vertices which make counting much more complex. 
Specifically, in a locally optimal set all vertices have the majority of their neighbors in their own set. Because 
of this complication a clean closed formula for the number of configurations that give a locally optimal set 
with c edges coming out probably does not exist. Given that, our main technical effort is concentrated on 
determining the degree distribution of maximum probability. Using a (different kind of) local optimality 
argument we show that the number of neighbors each vertex has on the other side of the partition can be 
assumed to follow roughly a binomial distribution, with some undetermined parameters (Lemma [2]). 

Having arrived at this technical tool we then apply it to the problem at hand. First, as a warm-up we 
show that if one ignores the degree constraints it is not hard to calculate the parameters of the degree 
distribution. We thus, after some calculations, arrive at exactly the same bound given in [S], albeit through 
a rockier road and gaining some valuable insight in the kind of configurations that lead to bad cuts. Then we 
go on to exploit this insight combined with the local improvements idea to improve on the expansion lower 
bounds for specific values of A. This leads to a set of optimization problems which we solve numerically to 
obtain our new results (note that the correctness of the bounds we give can however be verified by simple 
calculations). Our bounds improve on those given in [5] and for A > A improve the state of the art (for cubic 
graphs the best result is still given in [TT], but the complicated analysis there is restricted only to A = 3). 
Finally, we show how our analysis can be applied for general A and that an improved asymptotic lower 
bound on the expansion of random Z\-regular graphs can be obtained. The new coefficient of A we get offers 
only a small improvement over the one given in [5], so we do not expend much effort trying to calculate it 
exactly. Instead, for the sake of simplicity we simply show that a strictly smaller coefficient is achievable. 
Despite the lack of a large improvement in this constant, we believe there is still interest in this result, in 
that it indicates that the bound given by Bollobas is not the "right" bound, and also in that it shows that 
the local improvement method we employ here applies generally and not just for small values of A. 



2 Definitions and Notation 

We use standard graph terminology and will be dealing with regular graphs only. We use A to denote the 
degree, n to denote the number of vertices of a graph. When we talk about random regular graphs we mean a 
graph constructed via the following process. Take An vertices (assume An is even), numbered 0, . . . , An — 1. 
First, select a perfect matching on these vertices uniformly at random among all perfect matchings. Then, for 
all k G {0, . . . ,n} merge the vertices kA, kA + 1, . . . , (fc + 1)A — 1 into a single vertex of degree A. Observe 
that though this process may produce a multi-graph, the probability that it produces a simple graph is 
bounded away from as n — >■ oo and all simple Z\- regular graphs are equiprobable ([5i). The probability 
space therefore consists of {An — l){An — 3) . . . matchings. We denote this product by (Z\ri)!! (and generally 
we denote by n!! the product of all odd positive integers which are less than or equal to n). We use logn to 
denote the binary logarithm of n and ln?i to denote the natural logarithm of n. 
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Let S C V he a, set of vertices with \S\ < n/2. Let c{S) denote the number of edges with exactly 
one endpoint in S. The edge expansion of S is defined as ^j§j-- T^^^ edge expansion of G is defined as 
i{G) = min|s|<„/2 

Given a partition of V into S,V \ S we define the out-degree of a vertex as the number of neighbors the 
vertex has on the other side of the partition. We will characterize cuts mainly by the number of vertices of each 
out-degree in S and V\S. For this, we will use {A + l)-dimensional vectors, denoted as s = {so, si, . . . , sa)- 

More specifically, consider a set 5 C V, with \S\ < n/2. Let Si,i G {0, . . . , Z\} denote the number of 
vertices in S that have exactly i neighbors in V \ S. We write s — (sq, . . . , s^) and we have Si — \S\. 
Note that by definition — c{S). Also, we denote by s\,i £ {0, . . . , A] the number of vertices in \ 5' 

that have exactly i neighbors in S and let s' — (sq, . . . , s^). Note that is[ — c{S) and J^i + 12i — 
We use d to denote the maximum i such that Si > and d' to denote the maximum i such that > 0. 



3 Main Analysis 

In this section we establish the main tools we will need in the analysis of the expansion of random regu- 
lar graphs. The main result gives a general form for the distribution of out-degrees s that has maximum 
probability, conditioned on the size of the cut (Lemma [2]). In the end of this section we apply this tool to 
re-prove the main theorem of [5], giving the additional insight that the distribution of out-degrees in the 
configuration of maximum probability is a binomial distribution. In the following sections we use Lemma 
[2] together with some observations on d, d' to obtain improved results. However, for now we will work only 
with the basic case where there is no constraint on the out-degress and d ^ d' ^ A. As a result all sums 
given in this section can be read to range from i = to Z\. 

We are given a random Z\-regular graph G{V,E). We want to examine the asymptotic (in n) behavior 
of the expansion of such Z\-regular graphs. In other words, for each fixed A we would like to have a lower 
bound on a random graph's expansion which holds with high probability for sufficiently large graphs. For 
concreteness, say we always assume that n > A^'^ . To avoid inessential difficulties let us also assume that n 
is even. 

Fix a set of vertices S. Our first step is to write down exactly the probability that a certain configuration 
of out-degrees occurs for S and V \ S and that it produces a cut of size c{S) — c. More precisely, we will 
calculate the probability that, given s = (sq, si, . . . , S/i) and s' — (sq, si, . . . , s/i) such that — \^\^ 

^ . s'^ = \V\ S\ and isi = is- = c the set S has exactly c edges coming out, and these are distributed 
in a such a way that for all i there are Si vertices in 5' with i neighbors in \ 5 and s[ vertices in y \ S" with 
i neighbors in S. That probability is 

\S\\ -^fAV' \V\S\l VtM\'' ,{A\S\-c)ll{A\V\S\-c)U 



C! 



Let us explain this. The first factor comes from deciding how to partition the set S into groups of vertices 
with the prescribed number of neighbors in \ 5*. Recall that each vertex in our graph is derived from a 
group of A vertices in the original configuration where we take a random perfect matching. Therefore, if we 
decide that a vertex has i edges to vertices in y \ 5 we must decide which of the A vertices of the group 
are the endpoints of these edges. This gives the second factor. The third and fourth factors are similar for 

V \ S. The c! comes from deciding how to match the vertices from each side that have neighbors on the 
other. Finally, we take the product of the total number of matchings for the rest of the vertices of S and 

V \ S. This is divided by the total number of matchings overall. Recall now that we have written P{s, s') 
for a specific set of vertices S. By union bound, the probability that some "bad" set of vertices exists, such 
that its size is u, the cut has c edges and the out-degree distributions of the set and the complement are 
described by s,s' is at most (")P(s,s'). 

Let lLJ(u, c) be the set of all pairs (s, s') such that J2i — J^ii^i + ^i) — ''^ ^"^^ X^i — X^i ^^'i — 
In other words, U is the set of all possible vectors describing a configuration where liSj = u and there are c 
edges crossing the cut. Suppose we could show that 
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ri/2 L0"J / \ 

li-E E E nPis,s')^0 (1) 



u=0 c=0 (s,s')eO(«,c) 



for some parameter (j) depending on A. Then we would have that with high probabihty a random A- 
regular graph has expansion at least (j) from a simple application of the union bound. Our goal is to determine 
the largest possible parameter </)(Z\) such that this holds, and recall that by the results of [S] and [T] we expect 
the correct answer to be ^ ~ 0{^/A). 

Observe that the above sum has at most Z\n^'^+^ terms, since u < n, c < Au and there are at most n"^ 
different vectors s such that \s\ = u for all u. Let Sjvf,s^ be such that P{smtS'j^) is maximum among 
all {s^s') E U(u, c) for all u < n/2 and c < (j)u. It is sufficient for us to have that (2)-P(sAf, s^) = 

o{A^^n^^'^^^). In fact, we will have the much stronger condition {^)P{sm, Sjv^) — for some a > 1, 

or in other words we will show that the probability that a bad set exists is exponentially small. 

We have therefore reduced the question to that of determining a bad configuration with maximum 
probability and upper-bounding that probability. We would like to look for the vectors s, s', while avoiding 
some extreme cases. We can do that by using the fact that small perturbations in these vectors do not affect 
the probability too much. This is shown in the following lemma which essentially shows that changing some 
of the coordinates of a configuration vector (s, s') has a small impact on the probability. 

Lemma 1. If the configuration vectors {s, s'), {p,p') satisfy ma.Xi{\si — pi\} < ni and ma,Xi{\s^ — p^\} < ni 
we have P{s,s') < 2°(")F(p, p')- 

Lemma [T] allows us to neglect some annoying cases, such as for example sets with too small size {o{ni)). 
We now only need to find the configuration of maximum probability among configurations where Si > 
for all i and make sure that the probability is 2~^('^"). Then, by Lemma [1] we will know that looking for 
the maximum through the space of all configuration vectors would have made no difference. 

Now, given fixed values of u, c let us try to find the configuration vector s that gives maximum probability. 
It is not hard to see that we are trying to maximize the following function: 



Sq\si\ . . .Szl! 



In other words, we want to determine the distribution of the out-degrees of u vertices so that '^isi = c 
and F(s) is maximum. Note that the problem for the vector s' is essentially identical. 

Lemma 2. Let he a configuration vector such that sf^ = u, isf^ = c, for all i, sf'^ > ni and 
F{s^) is maximum among all vectors that satisfy the previous conditions. Then \sf^ — Sq^ ^ ^ (^)l — 
o{ni) for all i,Q < i < A. 

Proof. We define for i > 2 to be the vector with p^ = p\ = p\ = pI-i = ~1 and all other coordinates 
equal to 0. Intuitively, is a small perturbation vector, which does not change either the size of S or the 
size of the cut. Since maximizes the function F, it must be the case that both s^^ -f p* and s*^ — p' give 
smaller values. So we should have ^^^''Z , \■^ > f and also ^^^^ — \^ > 1. 

-l-P ) ~ -F(s — P ) ~ 

ihe hrst mequality gives — — }a\ > 1- ihe second mequality gives jnr, . ) i \ > 1- 

Observe that "^^^y-^ — 0(Z\^). Also, J/^m — *^("^) because all coordinates of s*^ are at least ni . Thus, 
the first inequality gives gMji_^ > 1 - O(^). The second gives ^m'^m' > 1 - O(^). Together 
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these imply that l-n— 14 }'a\ ~ 1| = 0(—^). Since the right-hand side is very small, intuitively what we 

^1 ^i-i y ^ j V 

wiU do is to assume it is essentially and solve the recurrence relation. 

More precisely, let p be the vector defined as follows: po ~ Sq^pi = sf^ and pi — ^^p'~^ ^^^'i for all 

i>2. It is then not hard to see (or verify by induction) that pi — Sq^ (^ty) (a)- 
We have that 



From this we get that (1-0 (^^))p2 < < (1 + O (^■^^)p2. By induction we then have (1 - 
O Yp, < sf < (1 + O Yp,. Using the fact that A < we have ((1 ± O = 1 ± O . 

Therefore, for all i we have |sf^ - Pi\ < Pi ■ O = 0{A^^/n) 0(713). □ 



Informally, as a result of Lemma [2] we can now assume that the vector describing the configuration of 
maximum probability has Si = Sq ^37^^ i^)' since for sufficiently large n we know that the maximum has 
an edit distance from this vector smaller than the one allowed by Lemma [1] To ease notation slightly we 
introduce the parameters /3,7 defined so that sq = /^X^i -^i ^^'^ ^1 ~ I'^so- We will also use /3',7' similarly 
as parameters corresponding to the vector s' . 

Going back to equation ([T]) we would like to show that the maximum probability is achieved for bisections, 
that is, u = n/2. This is intuitively unsurprising, and it is stated by the following lemma. 

Lemma 3. Let b,b' be two vectors such that "^^bi = X]j ~ '^/2; X^i*^* ~ Xi ~ '^'^'^ these two 
vectors maximize the quantity (T)P{b,b') among all vectors that satisfy the previous two conditions. Then 

for all vectors s,s' with J2^s^=u< n/2, Y^i^Si = = c we have (")P(s,s') < 2°(") (2)P(b, b')- 

Thanks to Lemma [3] we are only interested in partitions of size exactly n/2, since these give the maxi- 
mum probability of achieving a cut of size c (perhaps modulo a 2°^"^ factor which will prove insignificant). 
Therefore, for any expansion factor the probability of finding a set with expansion less than (j) is maximized 
when that set has size exactly ^. 

It is not hard to see that if the partition of the vertices is a bisection then the distribution of out-degrees 
that gives maximum probability is the same for the two parts. In other words, if s, s' are the two vectors 
that give maximum probability for a cut of size c and Xi — Xi ^'i then it must be that s — s' . Using this, 
the probability that a bisection which cuts c edges exists is upper bounded by 

(n/2)! y-AMy^'c!(^-c)! 



n/2/ V so!si! • • -szi! / -'■-'-V*/ {An)\\ 

where we have used the fact that (n!!)^ < n\. Now we would like to find the range of values of c for which 
the above quantity is very small, that is 2~'^('^"). To ease notation, we write c = (1 — rj)-^, where > 0, 
and look for the correct range of rj. 

Using the facts that for large n we have (^)" < n! < poly{n), (^^"2) < 2" and n!! > {^)^ poly{n) 
we get that 



log Pb (s) < n + n log 



2e 

(1 + ri)An 



A 

2 ^5, log- 



(1 - ■q)An (1 - r])An 



1=0 



4e 



(l + ri)An An, An 
log log— +o(n) 



5 



As mentioned, assuming that Si = /57'('^)§ cannot affect the probability by more than 2°("), so with 
some calculations and using the fact that = f ^isi = (1 — ^7)^ the upper bound becomes 

An 

log Pb <n — nlog/3 — (1 — ^)~2~ ~ 

+ - f log(l + r;) + ^ ^ log(l-77) + oH (2) 

What remains now is to find the values of /3,7 that give the maximum probability. Thankfully, because 
d ~ dl — A we can use the binomial theorem to find a clean solution to this problem as follows. Recall 
that E.s, = E./37^(f)t = + 7)^t - t- Also, E.^^. = E.*/37'(f)t = W(fpt = /37^(1 + 
^^^A-ln — _ Tj)^. Dividing these two equations gives 7 = and then /? — (^3^) ■ Plugging these 
values into the upper bound for Pb we get 



1^ < 1 - (1 - r?) log(l - r;) - (1 + ,7) log(l + r;) + o{l/A) 

The probability is therefore 2"^*^") whenever the right hand side is less than some arbitrarily small 
negative constant. 

Theorem 1. For each A, ifrj > satisfies (1 — 77) log(l — ?;) + (1 + log(l + 77) > then, for all sufficiently 
large n, random n-vertex A-regular graphs have expansion at least {1 ~ r/)^ with high probability. 

Thus, we have obtained again the main result of [5^, albeit through a much more winding road. One bonus 
is that through Lemma [2] we now have a better understanding of the out-degree distribution of maximum 
probability. Informally, one way to interpret this is to observe that by Lemma [2] the fraction of vertices 
that have out-degree i is proportional to the probability that a binomial random variable with parameters 
(Z\, takes the value i. This intuition will come in handy in the asymptotic analysis. 



4 Local Improvements 

In this section we introduce the idea of local improvements in the analysis. As mentioned, the main idea is to 
bound the probability that we find a partition that not only has small expansion, but also is locally optimal, 
in that flipping a pair of vertices cannot make the expansion smaller. The consequence of this property that 
we will rely on is that for a locally optimal set all vertices have out-degress at most (roughly) ^ ■ The analysis 
of the previous section will come in handy now, since the only difference is that d and d' are now smaller. 
Essentially all of it goes through unchanged, up to the point where we used the binomial theorem to find 
the values of /?, 7. Unfortunately, finding a clean formula for the values of /3, 7 in our case now becomes a 
much harder task. In this section we will be dealing with specific values of A, and we will therefore be able 
to calculate /3,7 and the minimum rj for which the graph is whp an (1 — 77)^ expander numerically. In the 
next section we will discuss ways to bound the values of /3, 7 to obtain an improved asymptotic bound as A 
tends to infinity. 

We wiU say that a set S' C j^l < n/2 is locally improvable if there exist u£S,v&V\S such that 

5 \ {u} U {v} has a strictly smaller cut. If S is not locally improvable we will say that it is locally optimal. 
Recall that d denotes the maximum i such that Si > and d' the maximum i such that > 0. Observe that 
whenever there exists a set with expansion (f) there must also exist a locally optimal set with the same size 
and expansion < 0. 

Lemma 4. If S is locally optimal then d+d' < A+1. Therefore, if S is locally optimal then mm{d, d'} < [-j] . 
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Proof. Suppose for contradiction that d + d' > A + 2. Let u E S have d neighbors in F \ S' and v E V \ S 
have d' neighbors in S. We swap u and v. First, suppose that u,v are not connected. The size of the new 
cut is c — d — d' + A — d+A — d'~c + 2A — 2{d + d') < c contradicting the local optimality of S. If u, v 
are connected then the new cut has size c — (d — 1) — {d' — 1) + A — d + A — d' < c, again contradicting the 
local optimality of S. □ 

In fact, we can say slightly more than what Lemma|3] states. Notice that the thorny case is when the two 
vertices of maximum out-degrees happen to be connected. If there are u E S,v E V\S with out-degrees d, d' 
respectively such that u, v are not connected then with the same argument we can say that if S is locally 
optimal then d + d' < A. Now, if there are at least A + 1 different vertices in S with out-degree d, then we 
can always find such a pair u, v and conclude that for a locally optimal set d + d' < A. Furthermore, if there 
are at most A vertices of out-degree d in S, then since A < n^o — o(n3) we can use Lemma [T] to say that 
the probability is essentially unchanged if we assume that there are no vertices of out-degree d, therefore 
the maximum out-degree in S becomes d — 1. As a result, we will from now on assume that d + d' < A for 
locally optimal sets. 

As mentioned, the analysis of the previous section goes through: even if vertices in S have maximum 
out-degree d and vertices in y \ 5 have maximum out-degree d' the probability of a certain configuration is 
still robust to small perturbations (Lemma [T]) and the most likely degree distribution is still in of the form 
Si = f3j^(^)\S\ (Lemma[2]). Also, bisections are again the interesting case, so we may assume that \S\ — ^. 

Revisiting the calculation that brought us to inequality ^ we now get: 

< 1 _ i log/3 _ i log/3' - (1 - 7j)j log7 - (1 - r;)^ logy - A 

+ log(l + .) + log(l - V) + oil) (3) 

Furthermore, we have 

'A\ /, /A 



: -E^v : =1 (4) 

E^n'ffl^E^v'ffUV^^ (5) 



= 1 ^ ^ i=l ^ 



Recall that in the previous section where we had d = d' = A we were able to use the binomial theorem to 
simplify the above equations and solve for /3,7,/3',7' as functions of ri,A. Plugging the result back into ([2]) 
and asking that the right hand side is negative led to Theorem [TJ Here, it is not clear how to do something 
similar. 

However, if A,d,d' have some known fixed values, we can do the following: pick a candidate value for 
rj and solve (numerically) equations (|4l5p for /3,7,/3',7'. Plug the solutions into ([3]) and check if the right 
hand side is negative. In such a case we can conclude that whp a random Z\-regular graph will not have a 
bisection with maximum out-degrees d, d' and '-''""'p'^" edges crossing the cut. Given that in a locally optimal 
set d + d' < A and allowing the maximum out-degree to increase can only increase the probability, we are 
essentially looking for the minimum rj such that the right hand side of ([3]) is negative for all possible pairs 
d, d' with d + d' = A. 

Performing this procedure for various small values of A we get the results listed in Table [TJ Despite the 
fact that (as with the bound given in [5 ) we need to rely on some numerical polynomial equation solver 
to find these numbers, their correctness can be verified relatively easily. Specifically, these bounds can be 
verified by checking that Yl'i=o Pl^ {^) — 1' Y^t=i P^l^ (f ) — ^~T^^^ ^^e same conditions apply for /3', 7' and 
the logarithm of the probability (as given in ([3])) evaluates to a negative number. 
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various values of A. The last column contains the 
in [5]. Due to space constraints we only list the final 



Observe that the improvements we obtain, though quite modest (in the order of 10^^ for Z\ < 10 and 
10~^ for somewhat larger A) seem to grow with A. This is a fact explained also by the results of the next 
section. 



5 Asymptotic Bound 

In this section we will extend the local improvement analysis to general A and show an asymptotic (in A) 
lower bound on the expansion of random Zi-regular graphs. As mentioned, this is interesting in part because 
it confirms the findings of the previous section that the gains from this analysis increase with A. 

Beyond this however, determining a tight asymptotic bound on this quantity is a fundamental mathe- 
matical question. Recall that the lower bound given by BoUobas is f - \/Z%/ln2. It is very natural to ask 
where the Vln 2 comes from and whether it is the "right" constant here. One intuitive explanation might 
be the following: consider a bisection of the set of vertices. Each edge of the random graph has probability 
(essentially) ^ of crossing the cut, therefore the expected number of edges crossing the cut is /i = If the 
edges were independent then the size of the cut would follow a binomial distribution with = ^ . What is 
the maximum size c such that we can guarantee that a cut of this size has probability o(2^") of occurring? 
If we approximate the binomial distribution by a normal distribution we get that "^Tf^^"^ must be 

o(2^"). Solving this gives c = f (^ — ^fA\/h\ 2). In other words, the analysis of [Sj shows that the true 
probability that a fixed set of vertices has c edges coming out can be approximated very well by assuming 
that edges are independent, an interesting and natural result. 

Nevertheless, as mentioned the weakness in the analysis of [5] is the union bound, which would only 
be tight if the probability that two different sets of vertices are expanding were independent. In some way 
the local improvement idea plays on this non-independence and, though we make only very small progress 
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towards bridging the gap between the known lower and upper bound, our results establish in an indirect 
way that \/ln2 is not the right constant and that there is something significant lost by the union bound in 
the analysis of .Si. 

To simplify presentation, in this section we define a locally optimal partition as one where mm{d, d'} < 
[y1 ■ This is a slightly weaker definition than what we use in the previous section. Assume without loss of 
generality that d < d' . Following similar reasoning as before, if there exists a non-expanding set there must 
exist one with d < [-y] , and since allowing d, d' to be larger can only increase the probability we can assume 
that d = Y and d' = A. Informally, we are applying the local improvement argument only to one side of the 
partition to simplify things. In particular, now we don't need to check various combinations of d,d', or to 
prove that the maximum probability is achieved when d = d' (which would agree with the numerical results 
of the previous section but would likely be complicated to prove directly). Finally, as usual assume that we 
are dealing with bisections, since these are the most interesting case and we will assume that rj < 
since otherwise we already know from [5 that the probability of such a set tends to 0. 

Having set d' — A we can again use the binomial theorem and equations (|4l5p to get 7' — j^f^ and 

/?' = (^^)^- Inserting these into ^ gives 



i^<l-|-^log/3-(l-r?)^log7 + o(l) (6) 
n 2 2 4 

The question therefore now becomes to determine /3, 7 and the smallest possible r/ so that the right hand 
side of ([6]) is negative. Going back to the left hand side of equations ()4|5|) we observe that 



We also have that 



i=0 ^ " ^ 1=0 



where we used the fact that ("^j — 4- {^^Z\) ■ 

The main intuition we are going to use now is that the two sums can be written more cleanly as proba- 
bilities using the binomial distribution. Specifically, let P\ be defined as P\ — Pr[_B(Z\, ::jqr[) < rf], where by 
B{n,p) we denote a random variable that follows the binomial distribution with parameters Also, let 
P2 be defined as P2 = Pr[B{A — 1, :^^) < d — 1]. Informally, as A tends to infinity we expect Pi and P2 to 
become almost equal. 

Now, with the help of the above equations (|4l5p become 

/3(7 + l)^Pi = 1 (7) 
f3j{j + l)^-^P2^^^ (8) 

Again, we will divide these two to get an equation for 7. To do this we would like to establish a relationship 
between Pi and Pa- Let P3 = Pr[B(Z\, = d]. 

Lemma 5. Let Pi, P2, P3 as defined above. Then Pi = P2 + ^^P^. 

To ease notation slightly let 6 — = ^ . Dividing equations (|7l8p we have — . Solving 

this for 7 gives 7 = j^^^p^- From this we get (3 — f ^2-26^ ) PT' ^^^^SS^^S these into (O we get 
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<l + ^ log(l + l log A - (1 - V)j log(l - V) - (1 + V)j log(l + 7? - 20) + 0(1) 
After re-writing all logarithms to base e we get 

< ^ + 2 ln(l - ^^) + I InPi - (1 - r;) ln(l - ,,) - (1 + r;) ln(l + - 20) + o(i) 

With a little foresight assume that in the end we will set r/ sufficiently large to have r] > 9 + -^^^ ■ In this 

case we would have = ^jE^ < ^{1 — r] + 9). Therefore, < y ~ Thus, the expected value of 
B{A, differs by at least half a standard deviation from d. We can use this to bound 9 from above using 

y _ 1 

the normal distribution to approximate Pi , P3 we have 9 <^J < . 

Now, using the fact that both rj and 9 are 0(^/^) we can use the expansion ln(l — x) = —x — ^ — 0{x^) 
and the above becomes 

41nPB 41n2 2 , o , 1 , 

< —r- + -rlnPi + -ry2 + - 

nA ~ A A ^A' 

So, we essentially need to have that rj^ > + ^ In Pi + 0^. if the sum of the last two terms is negative 

the bound for rj is strictly smaller that We can see that this is the case as follows. 

Recall that = ^j-^g > ^-y^- We are interested in the case 77 < ^^^^ . So, the expected value of 



B{A, differs from d = ^ by at most ^/ln2^/A which is at most 2 standard deviations. Again using a 
normal approximation for Pi we get that Pi < 0.84 InPi < —0.17. Using the upper bound on we have 
given already we have that the sum of the two terms is at most — Thus, the lower bound for 77 is strictly 
smaller than the one given in [5 . 



Theorem 2. There exists an a < 2-\/ln2 such that whp random A-regular graphs have expansion at least 

(l 9!-\A 



6 Conclusions 

In this paper we gave improved lower bounds on the expansion of random Z\-regular graphs for small values of 
A. We also showed that the asymptotic lower bound given by BoUobas is not tight. Perhaps more importantly, 
we gave an alternative analysis of the result of [S; incorporating the distribution of out-degrees and showing 
a general form. We believe that Lemma [2] could prove to be useful in establishing other properties of random 
regular graphs as well. 

As to the problems of this paper, much remains to be done. First, we would like to further refine our 
analysis, especially the asymptotic part. In particular, it would be nice to prove that when the local im- 
provement argument is applied to both sides of the partition the maximum probability is achieved when 
d = d' . This would agree with our numerical findings for small A, but a general proof has so far been elusive. 
Furthermore, it would be interesting to see if the local improvement idea can be taken any further, especially 
for special cases of small values of A where there may be room for more combinatorial arguments. Finally, 
for a concrete open problem, it would be nice to determine if there exist 5-regular graphs with expansion 
at least 1. This seems unlikely for random graphs, since the bound we give is still only (almost) 0.75 while 
there exists an upper bound of only slightly more than 1. However, if such graphs exist they would prove 
very useful in the field of inapproximability of bounded occurrence CSPs. 

Acknowledgement: I am grateful to Johan Hastad who first pointed out to me the idea of looking for 
a locally optimal set of vertices, which is the basis of this paper. 
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A Omitted Proofs 



A.l Proof of Lemma [T] 

Proof. We want to upper bound Since |s, ~ p,] < we have f4 < 20("''°s"). There are 2 A 

3 

factors of this form, so their total contribution is 2'^'"^"^ log") — 2°^"^. There are also 2 A factors of the form 

i'fY' ^\ each of which contributes at most 2"^"'' . Thus, their total contribution is at most 2^^"" — 2°^"\ 

using that A < . Similarly, the change in the remaining factors is also upper-bounded by 2^ "* — 2°("\ 
since the number of edges crossing the cut cannot change by more than A'^n* . □ 

A. 2 Proof of Lemma [3] 

Proof. Consider a pair of vectors s, s' as described above. We can assume without loss of generality that 
Si = /37*(f)u and similarly for s[ = I3'i'{'^){n - u), since by Lemma [2] and Lemma [1] the actual maximum 
cannot differ from the maximum achieved by such vectors by a factor of more than 2°'-"^. 
First, assume that Si — u < n/2 < s[. Therefore, we have 

E/^7^(f)-<E/^V'(f)(n-.) 

From these two we get that /3u < f3'{n — u) or Sq < Sg. Consider now the configuration that we can 
obtain from s, s' by decreasing Sq by one and increasing sq by one, thus increasing u. We claim that this new 
configuration achieves a value of (")P(s, s') that is at least as high. Indeed, because (") simplifies with the 
factors \S\\ = u\ and \V \ S\\ — {n — u)l into nl, the only factors affected by this change are , . But since 

So < Sq this product does not decrease by the change and we have a configuration with larger u. Repeating 
this argument leads to u = n/2. □ 



A. 3 Proof of Lemma [5] 

Proof. From the definitions of Pi, P2 it is not hard to see that Pi = P2 + Pr[P(Z\ — 1, :^^) — d — 1] ■ 
Intuitively, we will get at most d successes in A trials iff we either get at most d — 1 in the first A — 1, or we 
get exactly d in the first A — 1 and fail in the last. 

But now Pr[B{A-l, -r.)=d-l]-^ = (^-i)^^ = ^(^)^_^, where we have used the fact 

that(V) = ^(^). 
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